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ABSTRACT 

The statistics of gravitational lensing can provide us with a very powerful 
probe of the mass distribution of matter in the universe. By comparing pre- 
dicted strong lensing probabilities with observations, we can test the mass dis- 
tribution of dark matter halos, in particular, the inner density slope. In this 
letter, unlike previous work that directly models the density profiles of dark mat- 
ter halos semi-analytically, we generalize the density profiles of dark matter halos 
from high-resolution N-body simulations by means of generalized Navarro- Frenk- 
White (GNFW) models of three populations with slopes, a, of about -1.5, -1.3 
and -1.1 for galaxies, groups and clusters, respectively. This approach is an al- 
ternative and independent way to examine the slopes of mass density profiles of 
halos. We present calculations of lensing probabilities using these GNFW pro- 
files for three populations in various spatially flat cosmological models with a 
cosmological constant A. It is shown that the compound model of density pro- 
files does not match well with the observed lensing probabilities derived from 
the Jodrell-Bank VLA Astrometric Survey data in combination with the Cosmic 
Lens All-Sky Survey data. Together with the previous work on lensing probabil- 
ity, our results suggest that a singular isothermal sphere mass model of less than 
about 10^^h~^MQ can predict strong lensing probabilities that are consistent with 
observations of small splitting angles. 

Subject headings: cosmology:observations — cosmology:theory — gravitational lensing 
dark matter — galaxies:clusters:general — galaxies: halos 



1. Introduction 

Mapping the mass distribution of matter in the universe has been a major challenge 
for modern observational cosmology. The only direct procedure to weigh matter in the uni- 
verse is measuring its deflection of light by gravity. The statistics of gravitational lensing 
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can provide us with a very powerful probe of the mass distribution of the Universe. By 
comparing predicted lensing probabihties with observations, we can examine the mass dis- 
tributions of dark matter halos, in particular, their inner density slopes. It is well known 
that the Jodrcll-Bank VLA Astrometric Survey (JVAS) and the Cosmic Lens All-Sky Survey 
(CLASS) (Browne & Myers 2000; Myers et al. 2003; Browne et al. 2003) have provided us 
with observations of strong lensing probabilities for small image separations ranging from 
0.3" to 3". Based on the Cold Dark Matter (CDM) model, which has become the standard 
theory of cosmic structure formation, the lensing probabilities strongly depend on the den- 
sity profiles of CDM halos. The lensing model is usually described by a singular isothermal 
sphere (SIS), the Navarro- Prenk- White (NFW) model (Navarro et al. 1996, 1997), or gen- 
erahzed NFW (GNFW) density profiles of dark halos (Zhao 1996). Li & Ostriker (2002) 
employed a semi-analytical approach to analyze the gravitational lensing of remote quasars 
by foreground dark halos and checked the plausibility of various lensing models. They found 
that no model can completely explain the current observations: the SIS models predict too 
many lenses with large splitting angles, while the NFW models predict too few small split- 
ting angles. They therefore further developed a two-population halo model for lensing: small 
mass halos with a steep inner density slope and large mass halos with a shallow inner density 
slope, concluding that a combination of SIS and NFW halo models can reproduce the current 
observations reasonably well. 

Motivated by the sensitivity of the dependence of the image separation distribution of 
lenses below l" on both the inner mass profile of galactic halos and the faint end slope of the 
mass and luminosity functions. Ma (2003) compared the traditional approach that models 
lenses as SISs and the Schechter luminosity function, with another method that invokes a 
certain halo mass profile and the Press-Schechter mass function. She found that dark matter 
halos cannot all be SISs, otherwise, the mass function becomes steeper than the luminosity 
function and leads to a relatively high lensing rate on smaller angular scales. 

Li & Ostriker (2003) further proposed a lensing model based on three populations of 
halos distinguished by halos mass: Population 1 corresponds to normal galaxies of mass 
lO^'^/i~^M0 < M < IO^'^/i^^Mq with centers dominated by baryonic matter and inner density 
profile slopes of a = 2 (SISs); Population 2 corresponds to groups or clusters of galaxies of 
mass M > 10^^h~^MQ with centers of dark matter and a = 1.3 (GNFW models); Population 
3 corresponds to dwarf galaxies or sub-galactic objects with mass M < 10^^h~^MQ, whose 
centers lack baryons and are also dominated by dark matter, and also described by GNFW 
models of slope a — 1.3. Their results showed that both LCDM and OCDM cosmological 
models are marginally consistent with current lensing observations, while the SCDM model 
can be ruled out. 
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Besides these semi-analytical approaches, we can directly employ the results of high- 
resolution simulations for the density profiles of dark matter halos in the calculation of lensing 
probabilities. Jing & Suto (2000) performed a series of high-resolution N-body simulations to 
examine the density profiles of dark matter halos. They found a clear systematic correlation 
between halo mass and the slope of the profile at 1% of the virial radius, the slope of 
which is ~-1.5, -1.3 and -1.1 for profiles of galaxies, groups and clusters, respectively. They 
concluded that dark matter density profiles, especially in their inner regions, do not follow 
universal forms such as the NFW model. Therefore, consideration of the dependence of the 
slopes of density profiles on the masses of halos in calculating lensing probabihties seems to 
be necessary and reasonable. In this letter, we generahze Jing Sz Suto (2000) 's cosmological 
simulation results for density profiles of dark matter halos using the three-population GNFW 
models with slopes of « oc -1.5,-1.3 and -1.1 for galaxies, groups and clusters, respectively. 
This is an alternative and independent way to examine the slopes of density profiles of CDM 
halos. 



2. General Formalism for Lensing 

2.1. Lensing Equation for GNFW models 

The GNFW density profile can be expressed in the form p(r) = Psr^/r"(r -t- r^)^"" 
(Zhao 1996) where < a < 3. We generahze the inner slope from the simulation of dark 
matter halos using the GNFW model with a varying a which is dependent of the masses of 
halos. We take a = 1.5 for M < M^, a = 1.3 for < M < Me2 and a = 1.1 for M > Mc2, 
respectively, where Md ~ lO^^M© corresponds to the coohng mass scale (Porciani & Madau 
2000; Kochanek & White 2001) and Afc2 ~ 10^'^ Mq. The mass of a dark halo within r2oo can 
be defined as M = An /J^"" pr^dr = A-n p^rl f {ci) , and r2oo is the radius of a sphere around a 
dark halo within which the average mass density is 200 times the critical mass density of the 
universe. The function /(ci) = j^^ x'^dx/x°'{l + a;)^~" and ci = r2oo/rs is the concentration 
parameter 

Ci(Mi5, Z) = Cnorm^— ^[lOMl5]-°■'^ (l) 

where z is halo redshift, Cnorm = 8 (Bullock et al. 2001) and M15 = M/(lO^^h-^M0) is 
dimensionless halo mass. So Ps and Tg can be related to mass M15 and redshift z by 

^2/ x200 c? 1.626 M^l^ , , , 



-4- 



where pco is the critical mass density of the universe today and the expansion rate E{z) — 
A/riin(l + z)^ + ^A- The lensing equation for the GNFW profile is given hy y = x—iJLsg{x)/x] 
^ — XTs and f] — yTgd^/d^ are the position vectors in the lens plane and the source plane 

respectively, g{x) = udu {u'^ + z^)~"^^ {u"^ + z^'f^'^ + 1 dz^ and 



where [i^ is a parameter on which the efficiency of producing multiple images is strongly 
dependent. Here Scr = 4^ ^a^a is the critical surface mass density, and d^ = df^df^g/dg. 
dg and rff are the angular diameter distances from the observer to the source and to the lens 
object respectively, while df^g is the same quantity but from the lens to the source object. 
The lensing equation for a GNFW density profile with different a and setting /^^ = 1 is 
plotted in Fig.l. The curves are symmetrical with respect to the origin. Multiple images can 
be formed when \y\ <ycr, where t/cr is the maximum value of y when x < or the minimum 
value for x > 0. Generally speaking, there exist three images for \y\ < ycr- We will just 
consider the outermost two images stretched by the splitting angle when more than two 
images are formed. 



2.2. Lensing Probability for GNFW models 

We can write the cross-section as a (M, z) T^yl^'^l "& ( — A^o) with A^ > A^o in the 
lens plane for multiple images produced by a GNFW lens at z. ■»? is a step function and the 
splitting angle A^ is given by A^ = Vg^x/dj^ ~ 2xors/df^ where xq is the positive root of 
the lensing equation y{x) = 0. The lensing probability with image separations larger than 
A9 is given by Schneider et al. (1992) 

P{>A0)^ ^^^:^dz n{M,z)a{M,z)dM, (4) 
Jo dz Jo 

where Dp{z) = c/Hq dz/{l + z)E{z) is the proper distance from the observer to the 
lens at redshift z. The physical number density n{M, z) of virialized dark halos of masses 
between M and M + dM is expressed as n{M,z) = n{M,z){l + z)^, where n{M,z) = 
Pof{M, z)/M is the comoving number density (Press & Schechter 1974). po is the mean mass 

density of the universe today, and f{M,z) = -y^ MA dlnM ^^Ph ^^^^^ A'^{M) = 
2^ /o°° P{k)W'^{krM)k'^dk is the present variance of the fiuctuations within a sphere contain- 
ing a mass M. The power spectrum of CDM density fiuctuations P{k) — AkT'^{k) is given 
by Eisenstein & Hu (1999), where A is the amplitude normahzed to erg = A(rM = 8/i~^Mpc). 
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Fig. 1. — The lensing equation for a GNFW model with q;=1.5, 1.3 and 1.1 respectively. 
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Since the redshift distribution of quasars in the JVAS/CLASS survey is still poorly 
known, we adopt the mean value of Zg = 1.27 estimated by Marlow et al. (2000). In this 
letter, we will use spatially fiat ACDM models characterized by the matter density parameter 
fiin, vacuum energy density parameter Q\ and Hubble constant h = 0.75, and will calculate 
lensing probabilities with image separations greater than according to GNFW halo profile 
models with varying a. In the definition of cross-section or lensing probability, we will just 
consider the criterion A^, and neglect another one Qj. (the brightness ratio of the outermost 
two images) which is the ratio of the corresponding absolute values of the magnifications 
(Schneider et al. 1992). In order to investigate the effect of central black holes or bulges 
on lensing probabihty, Chen (2003a,b) introduced Qr into the calculation for lensing cross- 
section. Due to the existence of central black holes or galactic bulges, ycr becomes extremely 
large when |a;| approaches zero. Thus y^r can be determined by the consideration of 
together with A^^. However for GNFW halo models in the absence of central black holes or 
galactic bulges, we can see in Fig.l that the lensing equation curves are so smooth that we 
do not need to define cross-section by qr- 

Our numerical results are shown in Fig. 2 together with the observational ones, which 
for 6 < A^ < 15 from JVAS/CLASS takes the form of an upper limit. Based on the high- 
resolution simulation results (Jing & Suto 2000), we adopt a compound GNFW model for 
dark halos with three populations as described above. For the cosmological models, we first 
choose the new result from the Wilkinson Microwave Anisotropy Probe (WMAP): = 0.27, 
(78 = 0.84 (Bennett et al. 2003; Spergel et al. 2003) and M^i = IO^^Mq, the numerical result 
of which is depicted by a thin sohd hue in Fig. 2. It is clear that the predicted lensing 
probabilities are much lower than the observed values. Previous works (Chen 2003a,b) 
showed that lensing probability is sensitive to o"8 and Md- More specifically, it increases 
with l^in, o's and M^i- In order to improve the prediction, we take l^m) o's and M^i to be 
the higher values of 0.4, 1.2 and 5 x IO^^Mq respectively. As shown by the dashed line in 
Fig. 2, the derived lensing probabilities are still inconsistent with the observations although 
they have increased greatly compared to the former model. Thus we take the upper limit of 
1.6 for (78, as used in literature to date without changing the other parameters. The results 
are shown with a doted line in Fig. 2. There is almost no improvement in the small angle 
separations compared with the case of = 1-2, apart from the increase of the predicted 
probabilities on large angles. On the other hand, the lensing probabilities for both of cases 
as — 1.2 and as = 1.6 lie above the observed upper limit for large angle separations. 
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M^i = 1.0x10^^h"^MQ, 0^ = 0.27, ag=0.84 

'i'q M^i=5.0x10^^h"^MQ, 0^ = 040, ffg=1.20 

; M^^=5.0x10^-'h"^MQ, 0^=040, ag=1.60 ; 
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Fig. 2. — The lensing probabilities for image separations greater than using the compound 
GNFW dark halo profile models with varying a. The sohd-line histogram represents observed 
lensing probabilities from JVAS/CLASS (Chen 2003a; Chae et al. 2002). The thick sohd 
horizontal line denotes the upper limit of the null result for lenses with 6 " < < 15" 
from JVAS/CLASS. The other three curves are the predicted lensing probabilities in ACDM 
cosmology. 
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3. Conclusions and Discussion 

We have generalized the density profiles from high-resolution simulations of dark mat- 
ter halos according to GNFW models with varying a, and developed a compound model 
incorporating three populations of matter halos. We choose three typical cosmological mod- 
els: Model A, the new result of WMAP for model parameters: Qm = 0.27, erg = 0.84 and 
Mci = lO^^M©; Model B, in which = 0.4 is the upper limit adopted for the ACDM 
model, (78 = 1.2 and Md = 5 x lO^^M©; and Model C, which is the same as Model B 
except that as is taken to be 1.6, the maximum value used so far in the literature. We have 
compared our predicted results for lensing probabilities with JVAS/CLASS observations. 
Current observations for lensing probabilities from JVAS/CLASS only cover lensing events 
on small angle scales, though they also give the upper limit on separation angles over the 
range 6" to 15". Our numerical results as shown in Fig.2 show that the predicted lensing 
probabilities of the halos of three populations in current cosmological models do not agree 
with the JVAS/CLASS observations. For Model A, the predicted probabihties of lensing 
events lie below the upper limit observed on large angle scales, but is lower than the ob- 
servational results by a factor of about one order of magnitude on small separation angles. 
Even when the upper limit value of ACDM cosmology are used for fl^a and as, such as in 
Models B and C, the predicted lensing probability on small angle scales still do not match 
the observations. In the calculation of cross section, we have not considered the uncertainties 
such as the distribution of source redshifts for quasars, the brightness ratio g^j the scatter 
effect of the concentration parameter ci and the variation in Mc2- As discussed above, con- 
sideration of the brightness ratio Qr does not change our conclusion. As for the scatter in 
the concentration parameter ci, Chen (2003a,b) have investigated its effect on the prediction 
of lensing probabilities. By averaging lensing probabilities with a log-normal distribution, 
they found that the scatter in ci only increases the probabilities on larger image separations 
but has little effect at smaller separations. As for Mc2, any variation only contributes to 
the predicted lensing probabilities on large angle scales. In addition, for the sources in the 
JVAS/CLASS survey, the redshift distribution is still poorly understood except for the mean 
redshift < Zg >=1.27 as estimated by Marlow et al. (2000). Thus, most of the previous work 
(Li & Ostriker 2002; Chen 2003a) does not consider the redshift distribution of the sources. 
The prediction of Dunlop & Peacock (1990) model and the CLASS lensing sub-sample red- 
shift measurements suggest that the redshift distribution for CLASS unlensed sources can 
be modelled by a Gaussian distribution with mean redshift < Zg >—1.27. However this can 
not affect the conclusion of our letter. Even so, with the ever increasing number of observed 
lensing events, we look forward to the day when the redshift distribution derived from lensing 
observations can be taken into account for other aspects of lensing study. 

On the other hand, Klypin et al. (2001) have presented a convergence study of the 
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density profiles of cold dark matter halos simulated with varying mass and force resolutions. 
They concluded that, on radii larger than the "effective" spatial resolution, the density 
profiles of dark matter still do not experience any systematical trends even with the further 
increase of the number of particles or the force resolution. In addition, they also explained 
the systematic correlation between inner profile slope and halo mass found by Jing & Suto 
(2000), and pointed out that it is a trend in dark matter halo concentration. Navarro (2003) 
report recent results of numerical simulations designed to study the inner slopes of density 
profiles for CDM halos. Their results indicate that the inner slopes of CDM halos are 
considerably shallower than the asymptotic value of -1.5 proposed by Moore et al. (1999). 
This point probably indirectly strengthens our argument that the density profiles given by 
pure numerical simulation for CDM halos are inconsistent with the observational statistics 
of strong lensing. Therefore, in order to recover the observed strong lensing probabilities 
from the inner density profiles of mass halos, we can employ the SIS profile with a mass of 
less than about lO^^/i^^M©. This is in agreement with previous work such as Li & Ostriker 
(2002) and references therein. 
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